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Probing the N = 32 shell closure below the magic proton number Z = 20: 

Mass measurements of the exotic isotopes 52,53j^ 
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The recently confirmed neutron-shell closure at N = 32 has been investigated for the first time 
below the magic proton number Z = 20 with mass measurements of the exotic isotopes 
latter being the shortest-lived nuclide investigated at the online mass spectrometer ISOLTRAP. The 
resulting two-neutron separation energies reveal a 3MeV shell gap at A = 32, slightly lower than for 
®^Ca, highlighting the doubly-magic nature of this nuclide. Skyrme-Hartree-Fock-Boguliubov and 
ab initio Gorkov-Green function calculations are challenged by the new measurements but reproduce 
qualitatively the observed shell effect. 

PACS numbers: 21.10.Dr,21.30.-x,21.60.De,21.60.Jz 


Since the introduction of the shell concept in nuclear 
physics, the evolution of the shell structure far away 
from stability is one of the main research efforts. 
Similarly to atomic electrons, protons and neutrons can 
be regarded in a simple picture as occupying nuclear 
orbitals, of varying energy and angular momentum. 
Nuclides show spherical nature and enhanced stability 
when these orbitals are filled leading to the so-called 
“magic” numbers of nucleons, however the standard 
shell picture has been shown to no longer be applicable 
to some exotic species. As magic numbers disappear 
going farther from stability, new shell closures arise in 
consequence, as in the case of the doubly magic oxygen 
isotope [1, 2]. Such surprising isospin-dependent 
migrations of the nuclear quantum states point to the 
need to improve the understanding of the nuclear force 
[3-5]. 

High excitation energies of the lowest 2+ excited state 
are a traditional clue for closed-shell nuclides, however 
only for even-even species. Two-neutron separation 
energies (S' 2 „) are another strong indication - even in 
the case of odd-Z nuclides (such as potassium) and are 
obtained from mass measurements [6, 7] as performed 
with mass spectrometers like the ISOLTRAP experiment 
at ISOLDE/CERN [8]. 

Beta-decay studies of ®^Ca at ISOLDE published in 
1985 [9] revealed the energy of what was thought to be 
the first 2+ excited state (confirmed in 2006 [10]) and 


hinted at an enhanced N = 32 shell effect. Recent mass 
measurements including ®®’®'^Ca revealed a 4MeV drop 
in S 2 n between ®^Ca and ®'^Ca yielding a correspondingly 
large empirical shell gap [11]. An important test of 
the strength of the N = 32 shell closure comes from 
neighboring elements where proton correlations become 
active. Eor the isotopic chains of the heavier even-.Z 
elements titanium {Z = 22), chromium (Z = 24) and 
iron {Z = 26), the £'(2+) values have a local peak at 
N = 32 although with an absolute magnitude lower 
than for ®^Ca [12, 13]. Their two-neutron separation 
energies, however, as well as the ones of the odd-Z 
neighbor elements, which have no £’(2“'') indicator, show 
no significant kink at A = 32 - in sharp contrast to the 
closure of the Z = 20 proton core in calcium, where 
recently also for ®'^Ca a large £(2®') value has been 
revealed [14]. Until now, the region below Z = 20 was 
unexplored as these extremely exotic nuclides could not 
be produced in sufficient quantity. 

In this Letter, we present the first measurement of the 
masses of 52,53j^^ investigation of the 

N = 32 shell closure towards the neutron dripline for 
£ < 20. The experimental S 2 n values are compared 
to the predictions of self-consistent ab initio Gorkov- 
Green function (GGF) theory. This recent development 
significantly extends the reach of calculations based on 
realistic interactions, including three-nucleon forces, in 
particular to open-shell (and odd-£) nuclei. In contrast. 


2 


reference 
ion source 


D 



ISOLDE beam 


RFQ cooler and buncher 



pulsed 

drift 

tube 


precision Penning trap 

t 

preparation Penning trap 


MR-TOF 
mass separator 


movable 
MCP detector 


D=oai 


FIG. 1. Sketch of the part of the ISOLTRAP setup used for the measurements reported here. For details, see text. 


we show that the predictions of mean-field calculations 
with phenomenological interactions are highly dependent 
on the parameterization of both the particle-particle 
and particle-hole parts of the energy functional. To this 
end, we performed in addition Hartree-Fock-Bogoliubov 
(HFB) calculations with the Skyrme interactions SLy4 
and SLy5 [15], the latter containing a phenomenological 
tensor term. 

The exotic potassium isotopes were studied at 
ISOLDE/CERN with ISOLTRAP’s multi-reflection 
time-of-flight mass separator (MR-ToF MS) [16, 17]. Fig¬ 
ure 1 shows the parts of the setup relevant for the present 
measurements: a radio-frequency quadrupole (RFQ) ion 
cooler and buncher [18], the MR-ToF MS, and a micro- 
channel plate (MCP) detector. 

The nuclides and were produced by impact of 
1.4-GeV protons on a uranium carbide target. The atoms 
were surface ionized, along with stable chromium isobars. 
All ions were accelerated to 30 keV, transported via the 
ISOLDE magnetic high-resolution mass separator (HRS) 
to the ISOLTRAP setup [8, 19], captured and cooled in 
the RFQ buncher, and ejected towards the MR-ToF MS. 
The energy of the ions was adapted by a pulsed drift 
tube prior to the injection into the MR-ToF MS, where 
in-trap-lift capturing was applied [20]. After several hun¬ 
dred reflections the ions were ejected towards the MCP 
detector. 

Figure 2 shows time-of-flight spectra of isobaric ions at 
mass numbers A = 52 and A = 53. The calcium peaks 
were identified by an element-selective laser-ionization 
scheme. For ToF signals at longer flight times there are 
no known isobaric atomic ions except 
pected flight times of possible molecular ions produced 
at ISOLDE, such as oxides, dioxides, sulfides, hydrides, 
and of doubly charged ions have been compared with 
the observed peaks. The only species that would come 
close are ^®’^^Si(^®0)+, that contain short-lived silicon 
isotopes. However, due to the slow chemical release from 
the target, the probability for short-lived silicon beam 
components is very low. More importantly, it can be ex¬ 
cluded as only in the case of A = 52 the flight time of 
^®Si(^®0)+ would fit the measured peak; at A = 53, the 
measured mass deviates by 1000(130) keV from the value 
of ^^Si(^®0)''' [21], which corresponds to a time-of-flight 
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FIG. 2. Time-of-flight spectra for nuclides at mass numbers 
A = 52 (top) and A = 53 (bottom). Vertical lines indicate 
the Gr"*" fit ranges. 

difference in the order of 40 ns. 

In four hours, 12000 counts were collected for and in 
12 hours, 2300 counts for Note that with a half-life 
of just 30 ms, 52 k is the shortest-lived nuclide ever in¬ 
vestigated with ISOLTRAP. For mass calibration, stable 
from the ISOLTRAP reference ion source, as well as 
52,53(2;], jQj]g Kom ISOLDE were used. With the reference- 
ion masses mi ^2 and their corresponding flight times ti^ 2 , 
the ion mass of interest m is determined from its time of 
flight t by [11]: 

- my^) + {my + TO2^^)/2, (1) 

where the experimental values are described by the single 
variable CtoF = {‘2t — — t 2 )/ [2(ti — t 2 )]- Due to the 

deviation from Gaussian shape of the intense chromium 
ion peaks as shown in Fig. 2, the mean flight times 
have been determined from a limited region around the 
peak maximum (as indicated by vertical lines) that cover 
about 95% of the measured events. Uncertainties that 
arise from different choices of this region have been in¬ 
vestigated and taken into account. 

The resulting mass excesses ME = M — Au are listed 
in Tab. I, where M is the atomic mass, A is the mass 
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TABLE I. Half-lives ri /2 [21], reference nuclides (Ref. nucL), 
experimental ToF variable CtoF, and resulting mass-excess 
{ME) values of the investigated isotopes The mass 

values of the reference nuclides are taken from [21]. 

Isotope Ti /2 Ref, nucl._ Ctof _ ME (keV) 

110(4) ms ®®K,“Cr 0.50295425(256) -17138(33) 
30(5) ms 0.50306656(799) -12298(112) 


number, and u is the atomic mass unit. We note that, 
in principle, ®^Ca, with its mass well known from recent 
Penning-trap measurements [11, 22], would be the sec¬ 
ond mass reference of choice (in addition to ^^Cr) for 
measuring rather than from the reference ion 
source as ^^Ca appears in the same spectrum with 
and ^^Cr. In general, such a procedure (used here for 
the first time) is preferred because the MR-ToF data of 
all three isobaric ions are taken simultaneously, exclud¬ 
ing drifts between the measurements. The mass values 
from both calibrations are in agreement. However, the 
isobaric calibration has a larger uncertainty due to sig¬ 
nificantly lower statistics of the ®^Ca as compared to the 
signal from the reference ion source. 

In Figure 3, the two-neutron separation energies 
S 2 n = B{z, N) - B{Z, N -2), where B{N, Z) is the 
binding energy of a nuclide with Z protons and fV neu¬ 
trons, of the isotopic chains from argon to scandium are 
presented, including the present data. The new Sm val¬ 
ues show a significant drop from [22] to about 
3MeV, compared to only about 1 MeV from to 
clearly illustrating the N = 32 shell effect for potas¬ 
sium. The large uncertainty of the S 2 n value of ®^Sc 
(271 keV) [21, 23] at fV = 32 leaves room for some ambi¬ 
guity as to whether the S 2 n values of calcium and scan¬ 
dium cross. Thus, new precision-mass measurements of 
the corresponding Sc isotopes would be desirable. 

In Fig. 4 the empirical two-neutron shell gap 
S 2 n{N, Z) — S 2 n{N -f 2, Z) is plotted as a function of the 
proton number for N = 28 and N = 32, including the 
new point below the magic proton number Z = 20. A 
local maximum of the N = 32 shell gap has thus been 
revealed for the proton-magic Z = 20 core. In addition, 
the plot shows a similar behavior for both neutron num¬ 
bers going towards the dripline with lower Z, where the 
shell gap is rising up to a maximum at Z = 20 and falling 
by about 1 MeV at Z = 19. 

As shown already in [11] for the calcium isotopes, 
Density-Functional-Theory (DFT) calculations repro¬ 
duce the average S 2 n trend across N = 28 and N = 32, 
but do not predict the sharp drops associated with the 
two shell closures. We extend this study to Z < 20 
and perform Skyrme-Hartree-Fock-Bogoliubov (Skyrme- 
HFB) [25] calculations using the HFODD code [26] and 
interactions of the Skyrme-Lyon family [15], in the spher¬ 
ical approximation. 
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FIG. 3. Two-neutron separation energies for the Z = 18 — 21 
isotopes from the atomic-mass evaluation 2012 (AME2012, 
open symbols) [21], recent measurements of ‘^®’'*®Ar by Meisel 
et al. [24], and the ISOLTRAP data (from Ref. [11] and this 
work). 
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FIG. 4. Empirical shell gaps for A = 28 and N = 32 from the 
atomic-mass evaluation 2012 (AME2012, open symbols) [21], 
and the ISOLTRAP data (from Ref. [11] and this work). 


The calculations use the SLy4 or SLy5 interactions for 
the particle-hole part and a volume delta-pairing interac¬ 
tion for the particle-particle part of the energy functional, 
the latter being defined exclusively by a strength param¬ 
eter Vq. With respect to the SLy4 interaction, SLy5 con¬ 
tains in addition a tensor term. 

The top panel of Figure 5 shows the experimental and 
computed HFB S 2 n values for the potassium and cal¬ 
cium isotopic chains. The theoretical S 2 n are computed 
for nuclei of even neutron number. Self-consistent quasi¬ 
particle blocking of the odd protons is performed for the 
potassium isotopes, using the procedure described in [27]. 
A strength of the pairing interaction of — 200MeVfm^ 
reproduces the very smooth S 2 n trend observed in [11]. 
It describes correctly the experimental values on aver- 
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age but underestimates the drop at the crossing of the 
magic neutron numbers. A reduction of the strength of 
the pairing interaction (solid lines) leads to a significant 
improvement of the description of the experimental S' 2 „ 
trend. The addition of the tensor term with the SLy5 in¬ 
teraction leads to a change in the wrong direction. How¬ 
ever, a recent work [28] has shown that the effect of the 
tensor term in mean-field calculations strongly depends 
on the way it is constrained to experimental data. 

In addition to the empirical HFB approach, it is now 
possible to perform calculations up to the medium mass 
region using ab initio methods (see, e.g.. Refs. [29-36]). 
Thus, new mass calculations have been performed in 
the ab initio Gorkov-Green function (GGF) framework 
[31, 37, 38] that allows for the study of open-shell nu¬ 
clei. This method is particularly suited for the present 
purpose due to the ease of calculating odd-even systems, 
which also makes it a unique tool to investigate neigh¬ 
boring isotopic chains. 

In our calculations the only input are two- and three- 
body interactions fitted to properties of systems with 
A = 2, 3 and 4, without any further adjustments of the 
parameters. GGF calculations have recently addressed 
the region around Z = 20 [31] and are extended here for 
the first time beyond N = 32 for potassium. 

The present calculations made use of two- and three- 
nucleon forces derived within chiral effective field theory 
at next-to-next- and next-to-next-to-next-to-leading or¬ 
der (N^LO and N^LO), respectively [39, 40], extended to 
the low-momentum scale A = 2.0 fm“^ by means of free- 
space similarity renormalization-group techniques. The 
many-body treatment is set by a second-order truncation 
in the GGF self-energy expansion [37]. Model spaces 
up to 14 harmonic oscillator shells were employed and 
three-body interactions were restricted to basis states 
with Asmax < 16. Infrared extrapolations of the calcu¬ 
lated ground state energies were subsequently performed 
following Ref. [41]. We note that, in the present case, 
this procedure is formally defective due to the different 
truncations of one- and three-body model spaces. Nev¬ 
ertheless, we find that the trend expected from Ref. [41] 
is qualitatively reproduced, although with larger extrap¬ 
olation uncertainties. This is in agreement with other 
calculations [35]. As an example, we obtain binding en¬ 
ergies of 439.52(0.71) MeV for and 443.31(0.85) MeV 
for This overbinding of about 0.7 MeV/A is a gen¬ 
eral feature of currently available chiral interactions and 
it is a constant effect throughout the whole isotopic chain 
that cancels in separation energies [31, 35, 36]. 

GGF results for S 2 n of 47.49,51,53^; and 48.50.52.54^;^ 

shown in the bottom panel of Fig. 5 and are all result¬ 
ing from the infrared extrapolation. Different sources 
of uncertainty affect the present theoretical results (see 
Refs. [31, 38] for a detailed discussion). In particular, 
this method breaks particle-number symmetry (like HFB 
theory) and generates the correct expectation values for 
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FIG. 5. Two-neutron separation energies for the isotopic 
chains of potassium (left axes) and calcium (right axes), note 
the different scales. Open symbols: data from [21], filled sym¬ 
bols: calcium data from [11] and new mass data from this 
work. Top: with S 2 n values from HFB calculations using the 
SLy4 (green lines) and the SLy5 (red lines) interaction, with 
volume-type delta pairing of strength Vo = —150 MeV fm® 
(solid lines) or Vo = — 200MeVfm® (dashed lines). Bottom: 
with S 2 n values obtained from ab initio Gorkov-Green func¬ 
tion theory (see text for details). 


the proton and neutron numbers only on average, with 
a finite variance. However, the associated errors are ex¬ 
pected to cancel with good accuracy for energy differ¬ 
ences (such as S 2 n)- The errors indicated are uniquely 
those originating from the extrapolation fit and range 
between 0.4 to 1.5 MeV with increasing mass number. 
In general, GGF calculations are in fair agreement with 
measured S' 2 n, with the mismatch at being on the or¬ 
der of the truncation error. The significant drop from 
to is qualitatively reproduced but overestimated by 
theory, which also leads to an overestimation of the em¬ 
pirical shell-gap for potassinm. In contrast to the TV = 28 
gap, which is quantitatively better reproduced, the over¬ 
estimation of the N = 32 gap emerges as a common 
feature of ab initio calculations in this mass region (see 
also [35]). The present measurements therefore consti¬ 
tute an important test for the on-going developments of 
chiral interactions and ab initio methods. 

In summary, direct mass measurements performed 
with ISOLTRAP’s multi-reflection time-of-flight mass 
separator show that the neutron-shell closure at V = 32 
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extends below the closed proton shell Z = 20. The em¬ 
pirical gap is about 1 MeV weaker than for calcium indi¬ 
cating the stabilizing effect of the Z = 20 proton core and 
highlighting the doubly-magic nature of ®^Ca. Mean-field 
calculations for K and Ca show that interactions with 
effective parameters, although leading to more easily ac¬ 
cessible calculations, depend on the way the parameters 
are fitted and are thus not always precisely enough con¬ 
strained to yield robust predictions of shell effects. In 
contrast, Gorkov-Green function calculations, extended 
to > 32 along the potassium chain for the first time, 
offer predictions that are essentially parameter free and 
without any fitting to this mass region. By providing 
novel details on shell evolution, the measurements repre¬ 
sent an important test and are instrumental for validat¬ 
ing ab initio theories and extending them to heavier and 
more exotic nuclei. 

We thank T. Duguet and A. Gottberg for stimu¬ 
lating discussions. This work was supported by the 
German Federal Ministry for Education and Research 
(BMBF) (contract no. 05P12HGGI1, 05P12HGFNE and 
05P090DCIA), the Max-Planck Society, the European 
Union seventh framework through ENSAR (contract no. 
262010), the French IN2P3, the ISOLDE Collaboration, 
the DFG through Grant No. SFB 634, by the Helmholtz 
Alliance Program, Contract No. HA216/EMMI, by 
the United Kingdom Science and Technology Facili¬ 
ties Council (STFC) under Grants ST/L005743/1 and 
ST/L005816/1 and by the Natural Sciences and Engi¬ 
neering Research Council of Canada (NSERC), Grant 
No. 401945-2011. 

GGF calculations were performed using HPC resources 
from GENCI-CCRT (Grants No. 2014-050707 and 2015- 
057392) and the DiRAC Data Analytic system at the 
University of Cambridge (BIS National E-infrastructure 
capital grant No. ST/K001590/1 and STFC grants No. 
ST/H008861/1, ST/H00887X/1, and ST/K00333X/1). 
SK acknowledges support from the Robert-Bosch Foun¬ 
dation. 


[1] R. Kanungo et ai, Phys. Rev. Lett. 102, 152501 (2009). 

[2] C. R. Hoffman et al., Physics Letters B 672, 17 (2009). 


[3] D. Warner, Nature 430, 517 (2004). 

[4] O. Sorlin and M.-G. Porquet, Prog. Part. Nucl. Phys. 61, 
602 (2008). 

[5] O. Sorlin and M.-G. Porquet, Phys. Scr. T152, 014003 
(2013). 

[6] K. Blaum et al, Phys. Scr. T152, 014017 (2013). 

[7] L. Schweikhard and G. Bollen (eds.), Int. J. Mass Spec- 
trom., (special issue) 251(2-3), 85 (2006). 

[8] M. Mukherjee et al, Eur. Phys. J. A 35, 1 (2008). 

[9] A. Huck et al., Phys. Rev. C 31, 2226 (1985). 

[10] F. Perrot et al, Phys. Rev. C 74, 014313 (2006). 

[11] F. Wienholtz et al., Nature 498, 346 (2013). 

[12] M. R. Bhat, Evaluated Nuclear Structure Data File 
(ENSDF) (Springer Berlin Heidelberg, 1992) pp. 817- 
821. 

[13] R. V. F. Janssens et al., Phys. Lett. B 546, 55 (2002). 

[14] D. Steppenbeck et al.. Nature 502, 207 (2013). 

[15] E. Chabanat et al., Nucl. Phys. A 635, 231 (1998). 

[16] R. N. Wolf et al., Nucl. Instr. Meth. A 686, 82 (2012). 

[17] R. N. Wolf et al., Int. J. Mass Spectrom. 349-350, 123 
(2013). 

[18] F. Herfurth et al., Nucl. Instr. and Meth. in Phys. Res. 
A 469, 254 (2001). 

[19] S. Kreim et al, Nucl. Instr. Meth. B 317, 492 (2013). 

[20] R. N. Wolf et al., Int. J. Mass Spectrom. 313, 8 (2012). 

[21] M. Wang et al., Ghinese Phys. C 36, 16032014 (2012). 

[22] A. T. Gallant et al., Phys. Rev. Lett. 109, 032506 (2012). 

[23] H. L. Crawford et al., Phys. Rev. C 82, 014311 (2010). 

[24] Z. Meisel et al., Phys. Rev. Lett. 114, 022501 (2015). 

[25] P. Ring and P. Schuck, The Nuclear Many-hody Problem, 
Texts and monographs in physics (Springer, 2000). 

[26] N. Schunck et al.. Comp. Phys. Commun. 183, 166 

( 2012 ). 

[27] J. Dobaczewski et al.. Comp. Phys. Commun. 180, 2361 
(2009). 

[28] M. Grasso, Phys. Rev. C 89, 034316 (2014). 

[29] H. Hergert et al., Phys. Rev. Lett. 110, 242501 (2013). 

[30] A. Cipollone et al., Phys. Rev. Lett. Ill, 062501 (2013). 

[31] V. Soma et al., Phys. Rev. C 89, 061301 (2014). 

[32] G. Hagen et al, Phys. Rev. Lett. 101, 092502 (2008). 

[33] C. Barbieri and M. Hjorth-Jensen, Phys. Rev. C 79, 
064313 (2009). 

[34] G. Hagen et al., Phys. Rev. Lett. 109, 032502 (2012). 

[35] H. Hergert et al., Phys. Rev. C 90, 041302 (2014). 

[36] S. Binder et al, Phys. Lett. B 736, 119 (2014). 

[37] V. Soma et al., Phys. Rev. C 84, 064317 (2011). 

[38] V. Soma et al., Phys. Rev. C 89, 024323 (2014). 

[39] D. R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 
(2003). 

[40] P. Navratil, Few-Body Syst. 41, 117 (2007). 

[41] S. N. More et al., Phys. Rev. C 87, 044326 (2013). 



